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1 Introduction 

There is a prejudice among some specialists of non linear partial differential 
equations and differential geometry: convex analysis is an elegant theory but 
too rigid to address some of the most interesting and challenging problems 
in their field. Convex analysis is mostly attached to elliptic and parabolic 
equations of variational origin, for which a suitable convex potential can be 
exhibited and shown to be minimized (either statically or dynamically). The 
Dirichlet principle for linear elliptic equation is archetypal. 
Hyperbolic PDEs, for example, seem to be inaccessible to convex analysis, 
since they are usually derived from variational principles that are definitely 
not convex. (However, convexity plays an important role in the so-called en- 
tropy conditions.) Also, elliptic systems with variational formulations (such 
as in elasticity theory) often involve structural conditions quite far from con- 
vexity (such as Hadamard's "rank one" conditions). (However, convexity can 
be often restored, for example through the concept of polyconvexity [Ba] , or 
by various kinds of "relaxation" methods |Yol IABM] .) The purpose of the 
present paper is to show few examples of nonlinear PDEs (mostly with strong 
geometric features) for which there is a hidden convex structure. This is not 
only a matter of curiosity. Once the convex structure is unrevealed, ro- 
bust existence and uniqueness results can be unexpectedly obtained for very 
general data. Of course, as usual, regularity issues are left over as a hard 
post-process, but, at least, existence and uniqueness results are obtained in 
a large framework. The paper will address: 
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1. THE MONGE-AMPERE EQUATION 

(solving the Minkowski problem and strongly related to the so-called 
optimal transport theory since the 1990 's) 



(describing the motion of inviscid and incompressible fluids, interpreted 
by Arnold as geodesic curves on infinite dimensional groups of volume 
preserving diffeomorphisms) 

3. THE MULTIDIMENSIONAL HYPERBOLIC SCALAR CONSERVA- 
TION LAWS (a simplified model for multidimensional systems of hy- 
perbolic conservation laws) 

4. THE BORN-INFELD SYSTEM 

(a non-linear electromagnetic model introduced in 1934, playing an 
important role in high energy Physics since the 1990 's) 

Finally, let us mention that we borrowed the expression "hidden con- 
vexity" from a lecture by L.C. Evans about various models where the same 
phenomena occur (such as growing sandpiles |AEW] and weak KAM theory). 

2 Monge- Ampere equation and optimal transportation 
maps 

Given two positive functions a and (3 of same integral over R d , we look for 
a convex solution $ of the the Monge- Ampere equation: 



This nonlinear PDE is usually related to the Minkowski problem, which 
amounts to find hypersurfaces of prescribed Gaussian curvature. 

2.1 A weak formulation 

Assuming a priori that x G R d — > V$(x) is a diffeomorphism (with a ja- 
cobian matrix D 2 Q(x) everywhere symmetric positive), we immediately see, 
using the change of variable y = V$(x), that (JTJ is equivalent to the following 
"weak formulation": 



2. THE EULER EQUATION 



f3(V$(x))det(D 2 $(x)) = a(x), x e R d . 



(1) 
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for all suitable test function / on R d . In the words of measure theory, this 
just means that (3{y)dy as a Borel measure on R d is the image of the measure 
a(x)dx by the map x — > V$(a;). Notice that such a weak formulation has 
nothing to do with the usual definition of weak solutions in the sense of 
distribution (that does not make sense for a fully non-linear equation such as 
(jTj)). It is also weaker than the concept of "viscosity solution", as discussed 
in [Cg. 

2.2 A convex variational principle 

Proposition 2.1 Let us consider all smooth convex functions $ on R d with 
a smooth Legendre-Fenchel transform 

**(?/)= sup x ■ y - V(x). (3) 

xeR d 

Then, in this family, a solution $ to the Monge-Ampere equation (QJ) is a 
minimizer of the convex functional 

J[\ft] = J ty(x)a(x)dx + J V*(y)/3(y)dy. (4) 

Proof 

For any suitable convex function we have: 

J[*]= f ^(x)a(x)dx + j 'v*(y)0(y)dy = J(V(x) + ^*(V^(x)))a(x)dx 
(since V$ transports a toward (3) 

> J x ■ V$(x)a(x)dx 
(by definition of the Legendre transform (j3J)) 

($(x) + $*(V$(x)))a(x)dx 



(indeed, in the definition of 3>*(y) = sup x ■ y — $(x), the supremum is achieved 
whenever y = V<J>(x), which implies $*(V$(x)) = x ■ V$(x) — <3>(x)) 



$(x)a(x)dx + J $*(y)[3(y)dy = J[$], 
which shows that, indeed, $ is a minimizer for 
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2.3 Existence and uniqueness for the weak Monge- 
Ampere problem 

Based on the previous observation, using the tools of convex analysis, one 
can solve the Monge-Ampere problem in its weak formulation, for a quite 
large class of data, with both existence and uniqueness of a solution: 

Theorem 2.2 Whenever a and (3 are nonnegative Lebesgue integrable func- 
tions on R d , with same integral, and bounded second order moments, 



there is a unique L 2 map T with convex potential T = V$ that solves the 
Monge-Ampere problem in its weak formulation (TJj), for all continuous func- 
tion f such that: \ f{x)\ < 1 + \x\ 2 . 

This map is called the optimal transport map between a(x)dx and f3(y)dy. 

By map with convex potential, we exactly means a Borel map T with 
the following property: there is a lsc convex function $ defined on R d , val- 
ued in ] — oo,+oo], such that, for a(x)dx almost everywhere x G R d , $ is 
differentiable at x and V$(x) = T(x). 

Comments 

The usual proof |Br3[ ISK[ \RR\ IBr5j is based on the duality method introduced 
by Kantorovich to solve the so-called Monge-Kantorovich problem, based 
on the concept of joint measure (or coupling measure) |Kaj . However a 
direct proof is possible, as observed by Gangbo [GaJ. This theorem can be 
seen as the starting point of the so-called "optimal transport theory" which 
has turned out to be a very important and active field of research in the 
recent years, with a lot of interactions between calculus of variations, convex 
analysis, differential geometry, PDEs, functional analysis and probability 
theory and several applications outside of mathematics (see [VI] for a review). 
A typical (and striking) application to the isoperimetric inequality is given 
in the appendix of the present paper. 

3 The Euler equations 

3.1 Geometric definition of the Euler equations 

The Euler equations were introduced in 1755 [Eu] to describe the motion of 
inviscid fluids. In the special case of an incompressible fluid moving inside 
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a bounded convex domain D in R , a natural configuration space is the set 
SDiff(D) of all orientation and volume preserving diffeomorphisms of D. 
Then, a solution of the Euler equations can be defined as a curve t — > g t 
along SDiff(D) subject to: 

^o ft - 1 + Vft = 0, (5) 

where p t is a time dependent scalar field defined on D (called the 'pressure 
field'). 

As shown by Arnold [AKl EMJ, these equations have a very simple geomet- 
ric interpretation. Indeed, g t is just a geodesic curve (with constant speed) 
along SDif f(D), with respect to the L 2 metric inherited from the Euclidean 
space L 2 (D, R d ), and — Vpt is the acceleration term, taking into account the 
curvature of SDif f(D). From this interpretation in terms of geodesies, we 
immediately deduce a variational principle for the Euler equations. However, 
this principle cannot be convex due to the non convexity of the configura- 
tion space. (Observe that SDif f(D) is contained in a sphere of the space 
L 2 (D,R d ) and cannot be convex, except in the trivial case d = 1 where it 
reduces to the identity map.) 

3.2 A concave maximization principle for the pressure 

Surprisingly enough, the pressure field obeys (at least on short time intervals) 
a concave maximization principle. More precisely, 

Theorem 3.1 Let (gt,Pt) a smooth solution to the Euler equations |3|) on a 
time interval [io^i] small enough so that 

(ti - t ) 2 D 2 p t (x) < vr 2 , Vx G D (6) 

(in the sense of symmetric matrices). Thenp t is a maximizer of the CONCAVE 
functional 

Q I I q t {x)dtdx+ f J q [g to (x), g tl (x)]dx, (7) 

Jto JD JD 

among all t dependent scalar field q t defined on D. Here 

J q {x,y] = inf f\-q t {z{t)) + l -^^)dt, (8) 

where the infimum is taken over all curves t — > z(t) £ D such that z(t ) = x G D , 
z(t\) = y G D, is defined for all pair of points (x, y) in D. 
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Proof 

The proof is very elementary and does not essentially differ from the one 
we used for the Monge-Ampere equation in the previous section (which is 
somewhat surprising since the Euler equations and the MA equation look 
quite different). The main difference is the smallness condition we need on 
the size of the time interval. Let us consider a time dependent scalar field qt 
defined on D. By definition of J q : 



D Jq\gto(x),g tl (x)]dx < J t ]J<2^~dt^ - qt{gt{x)))dtdx. 

Using a standard variational argument, we see that, under the smallness 
condition the Euler equation §5§ asserts that, for all x G D 



M9t (x),g tl (x)] = jf (-\-^-\ 2 -p t (g t (x)))dt. 
Integrating in x G D, we get: 

/ J P \9to( x )>9ti(x)]dx = \ (-\—^\ 2 -p t (g t (x)))dtdx. 
Jd Jt Jd I dt 

Since gt G SDif f(D) is volume preserving, we have: 



Finally, 



(qt(x) ~ q t (gt(x))dx = J (p t (x) - p t (g t (x))dx = 0. 



q t (x)dtdx+ / J q [g t0 (x),g tl (x)]dx 



t JD JD 

< ( I p t {x)dtdx+ f J p [g to (x),g tl (x)}dx 

Jto JD JD 

which shows that, indeed, (pt) is a maximizer. 



3.3 Global convex analysis of the Euler equations 

The maximization principle is the starting point for a global analysis of the 
Euler equations. Of course, there is no attempt here to solve the Cauchy 
problem in the large for d > 3, which is one of the most outstanding prob- 
lems in nonlinear PDEs theory. (This would more or less amount to prove 
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the geodesic completeness of SDiff(D).) We rather address the existence 
of minimizing geodesies between arbitrarily given points of the configura- 
tion space SDif f(D). This problem may have no classical solution, as 
shown by Shnirelman phi J. Combining various contributions by Shnirelman, 
Ambrosio-Figalli and the author [Br4| ISh2[ IBr6l lAFj , we get the global exis- 
tence and uniqueness result: 



Theorem 3.2 Let go and g\ be given volume preserving Borel maps of D 
(not necessarily diffeomorphisms) and to < t\. Then 

1) There is a unique t dependent pressure field p t , with zero mean on D, that 
solves (in a suitable weak sense) the maximization problem stated in Theorem 

m\ 

2) There is a sequence g™ valued in SDiff(D) such that 

d 2 a n 

o {glT 1 + Vft - 0, 

in the sense of distributions and g$ — > g Q , g™ — ► gi in L 2 . 

3) Any sequence of approximate minimizing geodesies (in a suitable 
sense) betwween g and g\ has the previous behaviour. 

4) The pressure field is well defined in the space L 2 (]to, ti[, BVi oc (D)). 



Of course, these results are not as straightforward as Theorem 13.11 and 
requires a lot of technicalities (generalized flows, etc.). However, they still 
rely on convex analysis which is very surprising in this infinite dimensional 
differential geometric setting. Notice that the uniqueness result is also sur- 
prising. Indeed, between two given points, minimizing geodesies are not 
necessarily unique (as can be easily checked). However the corresponding ac- 
celeration field — Vpt is unique! It is unlikely that such a property could be 
proven using classical differential geometric tools. It is probably an output 
of the hidden convex structure. Let us finally notice that the improved reg- 
ularity obtained by Ambrosio and Figalli [AFj (they show that p belongs to 
L 2 (]to, ti[, BVi oc (D)) instead of Vp a locally bounded measure, as previously 
obtained in [Br6] ) is just sufficient to give a full meaning to the maximization 
problem. (A different formulation, involving a kind of Kantorovich duality 
is used in |Br4[ IBr6j and requires less regularity.) 
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4 Convex formulation of multidimensional scalar con- 
servation laws 

4.1 Hyperbolic systems of conservation laws 

The general form of multidimensional nonlinear conservation laws is: 

d 

d t u t + J2d i {F i {u t )) = 0, 
i=i 

where u t (x) EV C R m is a time dependent vector- valued field defined on a 
d— dimensional manifold (say the flat torus T d = R d /Z d for simplicity) and 
each Fi : V C R m — > R m is a given nonlinear function. This general form 
includes systems of paramount importance in Mechanics and Physics, such 
as the gas dynamics and the Magnetohydrodynamics equations, for example. 
A simple necessary (and nearly sufficient) condition for the Cauchy problem 
to be well-posed for short times is the hyperbolicity condition which requires, 
for all £ G R d and all v G V the m x m real matric 

8=1 

to be diagonalizable with real eigenvalues. For many systems of physical 
origin, with a variational origin, there is an additional conservation law: 

d 

d t (U(u t ))+J2di(Gi(u t )) = 0, (9) 
i=i 

where U and Gi are scalar functions (depending on F). (This usually fol- 
lows from Noether's invariance theorem.) Whenever, U is a strictly convex 
function, the system automatically gets hyperbolic. For most hyperbolic sys- 
tems, solutions are expected to become discontinuous in finite time, even for 
smooth initial conditions. There is no theory available to solve the initial 
value problem in the large (see [Daj for a modern review), except in two 
extreme situations. First, for a single space variable (d — 1) and small initial 
conditions (in total variation), global existence and uniqueness of "entropy 
solutions" have been established through the celebrated results of J. Glimm 
(existence) and A. Bressan and collaborators (well posedness) |Gf j IBB] . (Note 
that some special systems can also be treated with the help of compensated 
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compactness methods |Taj . without restriction on the size of the initial con- 
ditions.) Next, in the multidimensional case, global existence and uniqueness 
of 'entropy solutions' have been obtained by Kruzhkov |Krj in the case of a 
single (scalar) conservation law (m — 1). 

Theorem 4.1 (Kruzhkov) 

Assume F to be Lipschitz continuous. Then, for all m G L 1 (T d ), there is a 
unique {u t ), in the space C° (R + , L 1 (T d )) with initial value u , such that: 

dtUt + V ■ {F{ut)) = 0, (10) 
is satisfied in the distributional sense and, 

for all Lipschitz convex function U defined on R, the "entropy" inequality 

d t (U(u t )) + V-(Z(u t ))<0, (11) 
holds true in the distributional sense, where 

Z{v) = f V F'(w)U'(w)dw. 
Jo 

In addition, for all pair of such "entropy" solutions (u,u), 

/ \u t (x) — u t (x)\ dx < \u s (x) — u s (x) \ dx, Vt > s > 0. (12) 

J T d Jrpd 

This result is often quoted as a typical example of maximal monotone 
operator theory in L l . (For the concept of maximal monotone operator, we 
refer to [Brzl IABM] .) The use of the non hilbertian space L 1 is crucial. 
Indeed (except in the trivial linear case F(v) = v), the entropy solutions do 
not depend on their initial values in a Lipschitz continuous way in any space 
IP except fot p = 1. This is due to the fact that, even for a smooth initial 
condition, the corresponding entropy solution u t may become discontinuous 
for some t > and, therefore, cannot belong to any Sobolev space W 1,p (T d ) 
for p > 1 . 

4.2 A purely convex formulation of multidimensional 
scalar conservation laws 

Clearly, convexity is already involved in Kruzhkov's formulation fllOflll) of 
scalar conservation laws, through the concept of "entropy inequality". How- 
ever, a deeper, hidden, convex structure can be exhibited, as observed re- 
cently by the author |Brl0j . As a matter of fact, the Kruzhkov entropy solu- 
tions can be fully recovered just by solving a rather straightforward convex 
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sudifferential inequality in the Hilbert space L 2 . For notational simplicity, 
we limit ourself to the case when the initial condition u$ is valued in the unit 
interval. 

Theorem 4.2 (YB, 2006, JErfflj ) 

Assume Uq{x) to be valued in [0, 1], for x G T d . Let Y (x, a) be any bounded 
function of x G T d and a G [0, 1], non decreasing in a, such that 

u Jx) = f 1 l{Y (x,a) < 0}da, (13) 
Jo 

for instance: Y (a) = a — Uo(x). Then, the unique Kruzhkov solution to [TO]) 
is given by 

ut(x)= f l{Y t (x,a) < 0}da, (14) 

Jo 

where Y t solves the convex subdifferential inequality in L 2 (T d x [0, 1]): 

e d t Y t + F'(a) ■ V x Y t + d V [Y t ], (15) 

where r][Y] = if d a Y > 0, and i][Y} = +oo otherwise. 

Observe that Y — > F'(a) ■ V X Y + defines a maximal monotone 

operator and generates a semi-group of contractions in L 2 (T d x [0, 1]) [Brz] . 

Sketch of proof 

Multidimensional scalar conservation laws enjoy a comparison principle (this 
is why they are so simple with respect to general systems of conservation 
laws). In other words, if a family of initial conditions Uo(x, y) is non decreas- 
ing with respect to a real parameter y, the corresponding Kruzhkov solutions 
u t (x, y) will satisfy the same property. This key observation enables us to use 
a kind of level set method, in the spirit of Sethian and Osher [OS} IQFj . and 
even more closely, in the spirit of the paper by Tsai, Giga and Osher |TGO] . 
Assume, for a while, that u t (x,y) is a priori smooth and strictly increasing 
with respect to y. Thus, we can write 

Ut(x,Y t (x, a)) = a, Y t (x,u t (x,y)) = y 

where Y t (x, a) is smooth and strictly increasing in a G [0, 1]. Then, a straight- 
forward calculation shows that Y must solve the simple linear equation 

d t Y t + F'(a)-V x Y t = (16) 
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(which has Y t (x, a) = Yq(x —tF'(a), a) as exact solution). Unfortunately, this 
linear equation is not able to preserve the monotonicity condition d a Y > 
in the large. Sub differential inequality (j!5p is, therefore, a natural substitute 
for it. The remarkable fact is that this rather straighforward modification 
exactly matches the Kruzhkov entropy inequalities. More precisely, as Y 
solves (IT5|) , then 



can be shown to be the right entropy solutions with initial conditions Uo(x, y). 
For more details, we refer to |BrlOj . 



Our approach is reminiscent of both the "kinetic method" and the "level set" 
method. The kinetic approach amounts to linearize the scalar conservation 
laws as (jT6l) by adding an extra variable (here a) . This idea (that has obvious 
roots in the kinetic theory of Maxwell and Boltzmann) was independently 
introduced for scalar conservation laws by Giga-Miyakawa and the author 
jBrlt IBr2| IGMj . Using this approach, Lions, Perthame and Tadmor [LPT] 
later introduced the so-called kinetic formulation of scalar conservation laws 
and, using the averaging lemma of Golse, Perthame and Sentis [[GPS], es- 
tablished the remarkable result that multidimensional scalar conservation 
laws enjoy a regularizing effect when they are genuinely nonlinear. (In other 
words, due to shock waves, entropy solutions automatically get a fractional 
amount of differentiability!). On the other side, the level set method by Os- 
her and Sethian [OS} IQFj describes functions according to their level sets 
(here Y(t,x,a) = y). This is a very general and powerful approach to all 
kinds of numerical and analytic issues in pure and applied mathematics. An 
application of the level set method to scalar conservation laws was made by 
Tsai, Giga and Osher in [TGOj and more or less amounts to introduce a 
viscous (parabolic) approximation of sub differential inequality ffT5j) . Finally, 
let us mention that some very special systems of conservation laws can be 
treated in a similar way |Br8j . 




Remark 
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5 The Born-Infeld system 

Using convential notations of classical electromagnetism, the Born-Infeld sys- 
tem reads: 

B x (D x B) + D 
d t B + V x - 1 ; = 0, V-5 = 0, 

!\ + D 2 + B 2 + (D x B) 2 



D x (D x B) - B 
d t D + V x - 1 ; = 0, V-Z> = 0, 

'1 + D 2 + fi 2 + (D x 5) 2 



This system is a nonlinear correction to the Maxwell equations, which can 
describe strings and branes in high energy Physics |Bol IBIl [BDLLl iGi] . Con- 
cerning the initial value problem, global smooth solutions have been proven 
to exist for small localized initial conditions by Chae and Huh [CH] (using 
Klainerman's null forms and following a related work by Lindblad [Li] )• The 
additional conservation law 

d t h + V • Q = 0, 
where 



h = s]l + D 2 + B 2 + (D x B) 2 , Q = D x B. 

provides an 'entropy function' h which is a convex function of the unknown 
(D, B) only in a neighborhood of (0, 0). However, h is clearly a convex func- 
tion of B, D, and B x D. Thus, there is a hope to restore convexity by 
considering B x D as an independent variable, which will be done subse- 
quently by "augmenting" the Born-Infeld system. 

5.1 The augmented Born-Infeld (ABI) system 

Using Noether's invariance theorem, we get from the BI system 4 additional 
('momentum-energy' ) conservation laws: 

m + v Q®Q-B®B-D®D 1 l%h + v , Q = 

h h 

We call augmented Born-Infeld system (ABI) the 10 x 10 system of equations 
made of the 6 original BI evolution equations 

,BxQ + D^ n DxQ-B, , , 

d t B + Vx ^ = d t D + V x f = 0, 18 

n h 
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with the differential constraints 

V-£ = 0, V •£> = (), (19) 

together with the 4 additional conservation laws ffTTl) but WITHOUT the 
algebraic constraints 

h = yj\ + D 2 + B 2 + (D x B) 2 , Q = D x B. (20) 

These algebraic constraints define a 6 manifold in the space (h, Q, D, B) G 
R 10 that we call "BI manifold". We have the following consistency result: 

Proposition 5.1 (Y.B., 2004 WH ) 

Smooth solutions of the ABI system [T7^18il9\) preserve the BI manifold [20]) . 
Therefore, any smooth solution of the original BI system can be seen as a 
smooth solution to the ABI system p?]llgpigj) with an initial condition valued 
on the BI manifold. 



5.2 First appearance of convexity in the ABI system 

Surprisingly enough, the 10 x 10 augmented ABI system (I17lll8|ll9p admits 
an extra conservation law: 

d t U + V • Z = 0, 

where 

u(h , Q ,n,B)^-±- D2 + B2 ^ 2 



h 

is convex (and Z is a rational function of h,Q, D, B). This leads to the 
GLOBAL hyperbolicity of the system. 

Notice that the ABI system looks like Magnetohydrodynamics equations and 
enjoys classical Galilean invariance: 

(t, x)^>(t,x + ut), (h, Q, D, B) ->(h,Q- hu, D, B), 

for any constant speed «GiJ 3 ! 

For a large class of nonlinear Maxwell equations, a similar extension can be 
down (with 9 equations instead of 10) as in [SeJj- It should be mentioned 
that a similar method was introduced earlier in the framework of nonlinear 
elastodynamics with polyconvex energy (see [Daj ). 
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5.3 Second appearance of convexity in the ABI system 

The 10 x 10 ABI (augmented Born-Infeld) system is linearly degenerate 
[DaJ and stable under weak-* convergence: weak limits of uniformly bounded 
sequences in L°° of smooth solutions depending on one space variable only 
are still solutions. (This can be proven by using the Murat-Tartar 'div-curl' 
lemma.) Thus, we may conjecture that the convex- hull of the BI manifold is 
a natural configuration space for the (extended) BI theory. (As a matter of 
fact, the differential constraints V • D = V • B — must be carefully taken 
into account, as pointed out to us by Felix Otto.) The convex hull has full 
dimension. More precisely, as shown by D. Serre [Se2j . the convexified BI 
manifold is just defined by the following inequality: 

h > yjl + D 2 + B 2 + Q 2 + 2\D x B - Q\. (21) 
Observe that, on this convexified BI manifold ff2T]) : 

1) The electromagnetic field (D, B) and the 'density and momentum' fields 
(h, Q) can be chosen independently of each other, as long as they satisfy 
the required inequality (1211) . Thus, in some sense, the ABI system describes 
a coupling between field and matter, original Born-Infeld model is purely 
electromagnetic. 

2) 'Matter' may exist without electromagnetic field: B = D = 0, which leads 
to the Chaplygin gas (a possible model for 'dark energy' or 'vacuum energy') 

d t Q + V ■ (^f^) = V(k dth + V • Q = 0, 
h h 

3) 'Moderate' Galilean transforms are allowed 

(t, x) -> (t, x + U t), (h, Q, D, B) -+(h,Q- hU, D, B), 

which is impossible on the original BI manifold (consistently with special 
relativity) but becomes possible under weak completion (see the related dis- 
cussion on 'subrelativistic' conditions in [Br9] ). 

6 Appendix: A proof of the isoperimetric inequality 
using an optimal transport map 

In this appendix, we describe a typical and striking application of optimal 
transport map methods. Let Q be a smooth bounded open set and B\ the 



14 



Hidden convexity in nonlinear PDEs 



unit ball in R d . The isoperimetric inequality reads (with obvious notations): 

\Q\i-i/d\ B a/d < h dn \ 
a 

Let V$ the optimal transportation map between 

a(x) = ^l{xeQ}, (3(y) = j^l{xeB 1 }. 

In such a situation, according to Caffarelli's regularity result |Caj . V$ is a 
diffeomorphism between Q and B\ (up to their boundaries) with C 2 internal 
regularity (which is not a trivial fact) and 

det(D 2 $(x)) = iGil. 

holds true in the classical sense. Then the proof (adaptated from Gromov) 
of the isoperimetric inequality is straightforward and sharp. Indeed, since 
V$ maps Q to the unit ball, we have: 

\dtt\ = I dcr(x) > I V$(x) • n(x)da(x) 

JdQ JdVt 

(denoting by da and n{x) respectively the Hausdorff measure and the unit 
normal along the boundary of Q) 

= / A§(x)dx 
Jn 

(using Green's formula) 

> d I (det(D 2 $(x)) 1/d dx 
Jn 

using that (detAf/ d < l/d Trace(A) for any nonnegative symmetric matrix 
A) 

= d\Q\ 1 - 1/d \B 1 \ 1/d 

since det(D 2 <S>(x)) = ^ , iGfi. 
So, the isoperimetric inequality 

\n\ 1 - 1/d \B 1 \ 1/d < hdn\ 
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follows, with equality only when f2 is a ball, as can be easily checked by 
tracing back the previous inequalities. Notice that Gromov's original proof 
does not require the map T to be optimal (it is enough that its jacobian 
matrix has positive eigenvalues). However, the optimal map plays a crucial 
role for various refinements of the isoperimetric inequality (in particular its 
quantitative versions by Figalli-Maggi-Pratelli [FMPj . for example). 
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